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Abstract 



We present a procedure for recovering the conformal factor of an anisotropic con- 
Ph ■ ductivity matrix in a known conformal class in a domain in M" with n > 2. The 

I method requires one internal measurement, together with a priori knowledge of the 

r~| ' conformal class (local orientation) of the conductivity matrix. This problem arises in 

■ the coupled-physics medical imaging modality of Current Density Impedance Imaging 
(CDII) and the assumptions on the data are suitable for measurements determinable 
from cross-property based couplings of the two imaging modalities CDII and Diffu- 
sion Tensor Imaging (DTI). We show that the corresponding electric potential is the 

J> . unique solution of a constrained minimization problem with respect to a weighted total 

I variation functional defined in terms of the physical data. Further, we show that the 

■ associated equipotential surfaces are area minimizing with respect to a Riemannian 
. metric obtained from the data. The results are also extended to allow the presence of 

fSJ i perfectly conducting and/or insulating inclusions. 



Keywords: Anisotropic, Hybrid Problems, Conductivity, Diffusion Tensor Imaging, Current 
Density Impedance Imaging 



X 

c3 ; 1 Introduction 

Biological tissues such as muscle or nerve fibres are known to be electrically anisotropic 
(see e.g. [351 EZ])- In this paper, we consider the problem of recovering an anisotropic 
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electric conductivity a of a body f2 from measurement of one current J in the interior. Such 
interior data can be obtained by Current Density Imaging (GDI), a method pioneered at the 
University of Toronto ([T71 EH]) that used a Magnetic Resonance Imager (MRI) in a novel 
way. We also rely on the MRI-based Diffusion Tensor Imaging (DTI) method to determine 
the conformal class of a, as in the new DT-CD-II method recently introduced and tested 
experimentally in ([23]). Thus, we assume that the matrix- valued conductivity function is 
of the form: 



with cro(x) known from DTI, and with the cross-property factor c{x) a scalar function to 
be determined. This assumption is motivated by a number of physical studies which have 
shown a linear relationship between the conductivity tensor and the diffusion tensor (see e.g. 
[HI 123] and further references therein) . 

We show that, in dimension n > 2, the cross-property factor c(x) can be determined 
from knowedge of the current J in f2 and of the corresponding prescribed voltage / on the 
boundary dVt. In fact, the only internal data we require is the scalar function 



(with denoting the inverse of the matrix (Tq). This turns out to be the appropriate 
extension of the corresponding earlier result for isotropic conductivities ([32]), where the 
interior data was the magnitude \ J\. 

The method we will be presenting is based on the minimization of a weighted total 
variation functional defined in terms of a(x)and cro(x). (See Theorem 12.31 for the precise 
statement). 

More generally, we will show that when VL contains perfectly conducting and/or insulating 
inclusions, then knowledge of a, ctq and / determines the location of these inclusions (except 
in exceptional CclSGS^ , clS well as the function c{x) (hence also the anisotropic conductivity 
cr) in their complement. 

1.1 Background and Motivation 

Mathematical work on non-invasive determination of internal conductivity has focused largely 
on the classical method of Electrical Impedance Tomography (EIT). There have been major 
advances in the understanding of this nonlinear inverse boundary value problem (see [40j 
for an excellent review; in particular, see [131 HBj for recent results on recovering anisotropic 
conductivities in a given conformal class for the special case of admissible manifolds). It 
has also been shown that the EIT problem is exponentially ill-posed, yielding images of low 
resolution (see [TB| I2F|). 

In a new class of inverse problems (that includes the one studied here) one seeks to 
overcome the limitations of the reconstructions obtainable from classical boundary measure- 
ments by using data that can be measured noninvasively in the interior of the object. These 
are known in the literature as hybrid problems (also as coupled physics, interior data or 



a{x) = c{x)(7q{x), 



(1) 



a = (ctq J • J)'i 



(2) 
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multi-wave problems), as they couple two imaging modalities to obtain internal measure- 
ments. For overviews of such methods see [HI I2D]- For imaging the electric conductivity, 
there are several approaches that combine aspects of KIT with MRI: MREIT, CDII, Elec- 
tric Properties Imaging (see [311 EH] for recent reviews) or with ultrasound measurements: 
Acousto-Electrical Tomography ([HI El [21]), Impedance- Acoustic Tomography ([H]). 

The starting point for the method presented here is the measurement of one applied 
current J{x) at all points x inside a region VL. We briefly recall the influential idea of [T71 EH] 
for obtaining such interior measurements using MRI. The current J induces a magnetic field 
B{x). The component of B parallel to the static field of the imager can be determined at any 
point inside Q from the corresponding change in the phase of the measured magnetization 
at that location. By performing rotations of the object and repeating the experiment with 
the same applied current, all three components of B can be recovered, and J{x) is then 
computed using Ampere's law: 

J(x) = —V X B{x) 

where /xq is the magnetic permeability (essentially constant in tissue). For our purposes, it 
is important to note that this Current Density Imaging (GDI) method works equally well in 
anisotropic media, as no knowledge of the conductivity is needed for the determination of 
the current density J(x). 

Inside the body being imaged the electric potential u{x) corresponding to the voltage 
f{x) on the boundary solves the following Dirichlet problem for the conductivity equation: 

V ■ aVu = 0, xencW (3) 
u\dn = f 

where a is the (generally tensorial) conductivity of the material. In the case of isotropic 
conductivities, (i.e. scalar a) considered in [311 |33l [321 [3ll [29] i^i the absence of insulating 
or perfectly conducting inclusions one can replace a in the above equation using Ohm's law 
|J| = (t|Vu| to obtain the quasilinear, degenerate, elliptic, variable coefficient 1-Laplacian 
equation: 

V-(|J|^) = 0, xen. (4) 

The above equation was first introduced, with the above derivation, in the article [19] , where 
the Neumann problem was considered and examples of non-existence and non-uniqueness 
were given to explain that additional data was needed for determining the conductivity. In 
the article [31] it was shown that equipotential surfaces, namely the level sets of u{x), are 

2 

minimal surfaces with respect to the conformal metric |J|"-i/„, with J„ the nx n identity 
matrix; this result was then used to treat the Cauchy problem for equation . The Dirichlet 
problem for equation (|lj) can also have infinitely many solutions (see [32]). A solution around 
this difficulty was found in [32j, where the approach via the partial differential equation (|1]) 
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was replaced by the study of the variational problem for which it is the Euler-Lagrange 
equation. It was shown that the solution of [3] is the unique minimizer for this problem. We 
recall these results in the following theorem. 

Theorem 1.1. (U^)Let Q C M", n >2 be a domain with a connected C^'" boundary, a > 0, 
and let fi denote Lebesgue measure on Vt . Let (/, \ J\) e C^'°'{dVL) x C°(i7) with \ J\ ^ 
fi-a.e. be associated with an unknown conductivity a G C"(0). Then 

Ua = argmin { / | J| | Vf |/i((ia;); ^ Ioq = f o.iT-d /^({Vw = 0}) = 0} 
exists and is unique. 

Furthermore, a = ^ L°°{yt) is the unique C^ipt) scalar conductivity associated to 
the pair (/, | J|). 

A generalization of the above result was later obtained in the article [29], where the 
isotropic conductivity was shown to be determined from knowledge of | J| on the complement 
of open regions on which cr may be zero (in the case of insulating inclusions) or infinite (for 
perfectly conducting inclusions). We refer the interested reader to details in [29] . 



1.2 Statement of Results 

In this article we will extend the imaging method described above to the case in which 
the conductivity is anisotropic and known to be of the form a{x) = c{x)aQ{x) where c{x) 
is an unknown scalar function and ctq G C^{fl, Mat{M.,n)) is a symmetric positive definite 
matrix-valued anisotropic term, assumed known. 

We shall first prove an anisotropic analogue to Theorem 11.11 as a prelude to the more 
general results accounting for inclusions and less restrictive function spaces. For this, we will 
need to precisely define the class of data that arises from physical measurements. 

Definition 1.2 (First notion of admissibility). A triplet {f, ao, a) G H^d^l)xL'^{il, Mat(R,n))x 
L^{Q) shall be said to be admissible if there exists a c{x) G L'^{Q) such that 



a = (cTg J ■ J) 2 , 

where the current 

J = —C(TqVu 

corresponds to the potential u G H^{Q), the weak solution to the BVP 

V ■ (ccToVu) = 0, X eVL 



u\ -f 
Ian ~ J- 

The first theorem we prove concerns the minimization of the following functional: 

J^[v] := / a((ToVt; ■ Vv)^dfi (6) 



n 
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where id denotes the Lebesgue measure on Q. We present the following uniqueness result in 
section |21 

Theorem 1.3. Let Q C be a domain with a connected C^'" boundary, a > 0, and let 
the triplet (/, cro,a) G C^'°'{d^l) x Mat(R, n)) x be an admissible triplet as in 

Definition with \ J\ > Lebesgue —a.e. in Q. Denote by a = co"o € C°(fi, Mat(]R, n)) the 
unknown generating conductivity for this triplet. Then the following minimization problem 

argmin {J^[v]; v G W^'^ n C(n), ix{Vv = 0} = 0} (7) 

'"\afi=f 

has a unique solution where the functional is as in 

Furthermore, the unique Mat(R, n)) conductivity generating the current density J 

while maintaining the boundary voltage f is given by a = c(x)(7o(x) with the conformal factor 
c determined from the formula 



Following this we establish, in the remainder of section O the geometrical result that 
equipotential sets u^^{\) '■= {x; u{x) = X} Dfl are in fact minimal surfaces with respect to 
a certain Riemannian metric on Q which is defined in terms of of cro(x) and a{x). 

After the above preliminary results, we prove the following main uniqueness result in 
section m which allows for inclusions and more general function spaces. The precise statement 
requires an appropriate notion of admissibility, formally defined in section H] which involves 
some technical extensions of the criteria in Definition 11.21 It also requires an extension of 
the functional (E]) appearing in Theorem 11.31 A full definition of the proper generalization is 
postponed until section |3]but it takes the form of a weighted total variation \Dv\tf,, where 
the weight is defined in terms of a and ctq and where |-Df is a weighted distributional 
gradient discussed in section |3l 

Our uniqueness result also requires certain natural assumptions on the regions of perfect 
and zero conductance Ooo and Oq, respectively, as is discussed in greater detail in that 
section. Further, the zero set of the function a, denoted S, are also assumed to have a 
certain topological feature discussed in section [31 

Theorem 1.4. Let Q C M", n > 2, be a domain with connected Lipschitz boundary and 
let if,ao,a) G H^/'^{dn) x \ (Ooo U Oq), Mat(M, n)) x L°°{Vt) be an admissible triplet 

generated by an unknown conductivity a in the sense of Definition \4-l\ Then the potential u 
is a minimizer of the problem 

u = argmin{ / \Dv\^ : v G BViod^ \ S) and v\on = /}, (8) 
Jq 

and ifu is another minimizer of the above problem, then u = u in VL\{\J\ = 0}. Consequently 

o = -—^-—,<j,eC^{n\Z) 
[aoVu ■ Vn)2 

is the unique conductivity generating the admissible data triplet (/, (To,a) and Z is an open 
set discussed later which accounts for the inclusions. 



5 



With this shown, we lastly prove that level sets of such solutions as above minimize the 
area functional 

^(S) = j adS 

over a suitable space of variations of subsets of S C fi. 

An outline of this paper is as follows. In section |2] we present the main ideas building 
up to Theorem 12. 3[ Beginning with a formal statement of the problem, we describe the 
data of interest in our first definition of admissibility. We then proceed to establish the 
uniqueness result under this assumption. This section avoids the more technical details that 
arise from consideration of singular, perfectly conductive, and/or insulating inclusions and 
also considers potentials satisfying a more restrictive amount of smoothness. These more 
restrictive assumptions help clarify the main ideas that will be useful later. In section [3] we 
will introduce some tools from geometric measure theory needed to define an appropriate 
weighted space of bounded variation for our minimization problem. Most of the results 
presented in this section originated in the article [2]. In section H] we formulate a more general 
notion of admissibility and present the proof of our main uniqueness result. Theorem 14.21 
We then prove, in section [5l that equipotential sets minimize an area functional defined in 
terms of current density measurements. Finally, some of the technical facts on existence 
and uniqueness of solutions to a limiting form of the conductivity equation, suitable in the 
presence of inclusions with zero or infinite conductivity, as well as an equivalent optimization 
result in this setting, are briefly presented in section [Hi Section [7] presents conclusions and 
acknowledgments. 

2 Anisotropic Current Density Impedance Imaging 

In this section we present a simplified exposition of the main results of this paper , in order 
to illustrate the basic ideas used in the argument and to motivate the more general results 
to be presented later. We also use this section to briefly introduce some of the key geometric 
measure-theoretic concepts we will need and expand upon later; some excellent references 
thereon may be found in [IH [HI [ISl ED] • 

2.1 Uniqueness in the Variational Problem for Inversion 

Assume that the conductivity a is of the form c{x)o-o{x) with c(x), {ao)ij{x) G C"(fi), a > 0, 
c{x) > and (Tq symmetric and positive-definite throughout Q. 
Throughout the paper we will be using the notation 

{^,vU ■■= i^oO-v, leUo := ((^oO -0^ e c^{n,w) (9) 

to denote the inner product induced by o"o, and the corresponding norm, where ■ will always 
be taken to denote the Euclidean dot product. We also define the space as 

Wl'\n) := {v G W^'\n), Vv = 0}) = 0} 
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In what follows V denotes the usual (i.e. non-covariant) partial differentiation and we use 
the Einstein summation convention over repeated upper /lower indices. 

We begin by showing that the solution u to the BVP dSj) is a minimizer of an action 
on Q that is defined in terms of the internal density magnitude This generalizes the 

corresponding result for isotropic conductivities in 



Lemma 2.1. Assume that {f,aQ,a) is an admissible triplet and let u be a solution to the 
corresponding forward problem Then u is a minimizer of the action integral ■ ] defined 
by the following 

J^[v] := / a{x)\Vv\^^dx, (10) 



i.e. the relation 

T[v] > T[u] (11) 
Holds for all v G VI/^'"'^(fi) satisfying v = f. 



Proof Let v G Wl'\Q). Since a comes from an admissible triplet, there is a choice of c( 



x) 



such that a(x) takes the form a = \J\ -i. Then 



Jn 

c{x)\Vu\^f,\Vv\„gd^ 
> I cix)\{Vu,Vv)^,\dii (12) 
aVu ■ Vvdjj, 

n 

fAs (13) 
an dn 

= - [ fJ-ndS 
Jan 

with n an outer-oriented normal to dQ and where, in line (fT3|l . we have integrated by parts 
and applied the conductivity equation on u. Equality holds in line ( lT2ll if and only if Vu 
and Vf are parallel fi — a.e. In particular, we have 



J^[u] = - fJ- ndS 
Jan 

which, on comparing with the above, shows that m is a minimizer, as claimed. 

□ 

In order to prove the main result of this section we shall need to recall some basic notions 
from geometric measure theory. Firstly, by 'H°'(S) we denote the d-dimensional Hausdorff 
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measure of a set S C f2 defined as 

n 

W^iT) := liminf{^(diamEjy, |J Ej D S, diamE^ < 5} 

3=1 jeN 

The super-level set of a non-negative function u{x) G ly^'^(fi), given by Et '.= Qr\{u > t} has 
so-called locally finite perimeter, in the sense that the vector-valued Radon measure VxEt 
satisfies IVxeJ < oo for almost all t. For such sets we shall be concerned with the reduced 
boundary. 

Definition 2.2. The reduced boundary d*E of a set with locally finite perimeter is the set 
of points in M" for which the following hold; 

1. For all e > one has /g^^^) IVxeI > 

2. The measure-theoretic outer normal determined by 

v(x) := — lim 

exists, and satisfies \i'{x)\ = 1. 

For a super-level set Et the unit normal utix) exists T-L^~^ — a.e x G d*Et (see the remarks 
in [32]). 

We now present the main result of this section. 

Theorem 2.3. Let Q C M", n > 2, be a domain with a connected C^'°' boundary, a > 0, 
and let {f,ao, G C'^'°'{dn) x CiU, Mat{R,n)) x C"(fi) be an admissible triple with 

\J\ > fi — a.e. in Q. Denote by a & the unknown generating conductivity for this 

triplet. 

Then the following minimization problem 

argmin{J^[v];v E Wl'^ f] CiU)} (14) 

■>'\dn=f 

has a unique solution 

Further, the unique C°(r2) conductivity generating the local current density J while main- 
taing a boundary voltage f is given by a = c(x)cro(x) with the conformal factor c determined 
by 

UL-i 



c{x) 



|V?v| 
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Proof. The proof is essentially the same as the proof of Theorem 11.11 given in [32] so we 
present a self-contained but tersely abbreviated proof and refer the reader to that article 
for full details. Note first that since the triplet is assumed admissible the previous lemma 
ensures that argmin^^^^^j{J^[u\; u G PV^jJ:'^ nC(f2)} is nonempty. We choose a minimizer and 
call it Ua- 

To show uniqueness, assume to the contrary that another minimizer to problem [TU say 
u G W^''^ n C{Q), exists. Recalling the proof of Lemma (2.11 one sees the Cauchy inequality 
used in ( 1T2|) ensures that Vmo- = X{x)Vu for some non-negative A /i — a.e.. From this it 
follows that 

Vmct _ Vm 
|Vmo-| IVm] 

holds /i — a.e. We show that this implies equality of the minimizers away from Lebesgue- 
negligible sets. 

In view of Lemma 2.2 of [32], the above yields an identification of measure-theoretic 
normal vectors a.e.. It follows that the super-level set Et = {u > t} H Q has a measure- 
theoretic normal Ut{x) = which is continuously extendible from the reduced boundary 
d*Etr\Q to the topological boundary dEt fl Q. It then follows from a result of De Giorgi (e.g. 
4.11 in [15]) that, for almost all t, the region dEt n f2 is a C^-hypersurface with unit normal 
z/((x). From this, a parameterization of each connected component of this hypersurface shows 
that u„ remains constant on each such set. 

We verify that, for each t such that dE^ is a C^-hypersurface, each connected component 
Uf of dEt intersects dQ. Indeed, if not, then the Alexander duality theorem ([27]) implies 
that Ht U dil admits a decomposition of the form 11^ U dVl = Oi U O2 U (R"/fi) with Oj 
open and connected. We claim that {dOi U 9O2) fl 11^ 7^ 0. Indeed, were this not the 
case then dOi fl dQ 7^ for i = 1,2 and we could again apply the Alexander duality 
theorem to obtain the decomposition of W^/Ut into a union of two bounded and unbounded 
open, connected regions, 11^ = O;, U 0„. Since, by the contrary assumption Oi U O2 U is 
connected and unbounded we have that Oi U O2 U ^7"^ C Ou and thus, on taking complements, 
Ob C (Oi UO2 Ul]'^)^ C Ut. This contradiction shows that {dOiUd02) nUt 7^ 0, as claimed. 
It thereby follows that 11^ fl dQ 7^ 0. In other words, connected components of almost all the 
level sets dEt reach the boundary dQ. 

Now define G := {t e M+ 

■ ^ \dEt — l^-Bt^ ^ ^+ "^^^ suppose there exists a ball 
B C Q whose closure is contained in Q and such that B H {u E G} = If so then since 
I Vttl 7^ fi — a.e., it must map B to the closed interval [a, b]. However fi{[a, 6]) > whereas 
/i(G'^ n fi) = 0, a contradiction. We conclude from this that the union of all preimages of G 
under u and u is dense in Q and therefore u = u ^ — a.e., as was to be demonstrated. 

Finally, with J = — cctoVmo- we have {ctq^J ■ J)^ = (c^ctoVmo- ■ V-Uo-)2 /x — a.e.. This gives 
the desired formula for c(x). 

□ 
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2.2 Equipotential Sets are Minimal Surfaces in a Riemannian met- 
ric Determined from the Data 

We close this section with some interesting geometrical results about the level sets of solutions 
to ([5]). Given (Tq and the magnitude | of the current, we define a Riemannian metric on Q 
and show that the level sets of the corresponding potential function have zero mean curvature 
in this metric. In section [5] we will prove the stronger statement that these equipotential sets 
are in fact area minimizing. These are generalizations to anisotropic conductivites of results 
proved in [331 EI] for the isotropic case. 

As is customary, we denote | A| := det A for A G Mat(M, n) (which should not be mistaken 
for the norm \V\„g of a vector field V, as we hope shall be clear from the context). 

Proposition 2.4. Let Q C M", n >2 be a domain with Lipschitz boundary andu G C^'°'{Vt), 
a > 0. Assume the conductivity o is of the form c{x)crQ{x) for c, ctq ^ C°'{Q) with ctq a known 
positive- definite matrix-valued function and that \ Vu\, c{x) > ^-a.e. where u is the potential 
corresponding to the conductivity a and current density J via J = —aVu. 
Define the following Riemannian metric Qij on Q: 



9ij ■■= iWo\\J\l-i)--'i(^o%- (15) 

Then inside Q one has that 

V . iVWljr^) - 0. 

\\9 Vn|L 



Proof. We begin by noticing that |(To| "^H^l "-i ^ ~ c^"*""-! "-i {|(j|(crVn ■ Vu)}"-^a ^ 
whereby, with the above choice of gij one has that 

g'^ = {|cr|((TVn ■ Vm)}i^(J 

n 

Defining m{x) := |cr|(crVM ■ Vu) it immediately follows that \g\ = ^^J^ . Since [[(^"^VMllg = 
{{g^^Vu) ■ g{g^^'Vu)}'^ we have ||(7""'^Vm||j, = \/ (5'~^Vm) ■ Vu. Then 

n+l ]_ 



'u 



V 



^Jm{x)(T'S/u.^ 



^Jm{x) 

It follows from the fact that u solves the conductivity equation that 

v.(vH/^)^o. 

\\9 Vu|L 



□ 
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The above result immediately implies the following. 

Corollary 2.5. Suppose that u,c,aQ are as is in proposition i2.4\) - Then the level sets ofu, 
u~^{X) := {u = X} r\Q are surfaces of zero mean curvature in the metric 

9^j = iWo\\J\l-i)^'K%. 



Proof. As in the preceding proof of Theorem 12.31 the level sets u~^{X) are C^-hypersurfaces 
for yU— a.e. A. The unit vector rj := ^^Jj^" is (yf-orthogonal to each such level set u~^{X) since 
if ^ G Tc-u"^(A) then ^ ■ Vm \^(zu-i-(x) = therefore 



|Vm| 
IVmI 
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g 







Thus ±?7 are unit normals to the hypersurfaces u~^{X) in the Qij metric. Since the mean 
curvature of a surface with unit normal n is H = divg{n), with divg being the metric diver- 
gence, we conclude from Proposition (12.41) that when u satisfies the conductivity equation, 
we have H = 0. 

□ 



3 Preliminaries for the General Case 

In this section we prepare to expand upon the results in the preceding section by considering 
the conductivity equation over domains which may contain insulating or perfectly conduct- 
ing inclusions, i.e. regions of infinite or zero conductivity, respectively. We shall give the 
appropriate reformulation of the forward problem (|5]) in this setting. We also discuss inte- 
gration by parts and coarea formulae for spaces of bounded weighted variation which will 
play a key role in our main general uniqueness result. 

3.1 Weighted Total Variation 

We start by presenting some needed preliminary results about functions of bounded weighted 
variation. We will always use the notation Xa{x) to denote the characteristic function of 
a set A. Often, we will abbreviate vectors and matrices in component form, in addition, 
as earlier, we will employ the Einstein summation convention of implied summation over 
repeated upper and lower indices wherever appropriate. 
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Let n C M"" be a bounded open set with connected Lipschitz boundary and let a G L°°{Q) 
be non-negative. While the function a is now allowed to vanish, we require that its zero set 
S" := {x G n : a{x) = 0} always be assumed to satisfy the following structural hypothesis 

5:= our, (16) 

where F is a set of measure zero with at most countably many connected components, 
'H'^~^{dQriS) = 0, and where O is a mutually disjoint union of finitely many C^-diffeomorphic 
images of the unit ball, possibly empty. 

We recall the space of functions of locally bounded variation on a subset S C is given 

by 

BVioci^) '•- ^ LipilT), I \Du\ < oo, C E, K compact} 

Jk 

where Du is the distributional gradient of u. This generalizes the space BV{VL), the space 
of all L^iVL) functions with bounded total variation of the distributional gradient, i.e. those 
functions satisfying 

\Du\ < oo 



/ 



Let ctq G C^{VL^ Mat{W^n)) be a symmetric positive definite matrix with components 
((To)ij satisfying 

n 

< $^(ao).,(x)ee < Mier vx g ^ \ 5, ve g m", 

for constants < m, M < oo. We then denote by ^) the following function 

n 

yp(x,0 = a(x)(5^(ao),,fe)^. (17) 

For u G BVioci^ \ ^) define the weighted total variation of u, with respect to ip, in Q as 

\Du\^ = sup / uV ■ B (i/i, (18) 
Be<Ba Jn 

where 

= {B e L^{n,W) : V-fiGL"(fi) and \B\^-i < a{x) a.e. in fi}. 

and L^(fi,]R") is the space of vector fields of compact support in fl whose components are 
in L°^ln). 

n 

It is a straightforward consequence of the definition (fTSi) that \Du\ip is L['^~^^ (i7)— lower 
semi-continuous. It was shown in |2] by Amar and Bellettini that for any u G BV{Q), one 
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has the following integral representation formula for the weighted total variation appearing 
in equation (ITS!) . 

/ \Du\^ = I h{x,v'')\Du\ (19) 
Jn 

where, in the above, 



a Jn 



hix^v"") := {\Du\- ess sup {B ■v'')){x) \Du\ - a.e. x e n, (20) 

BG5B 

and denotes the vectorial Radon- Nikodym derivative v^{x) = j^^- Note that the right- 
hand side of equation ( !T9|) makes sense, as is iD-ul-measurable, and hence h{x,v^{x)) is 
as well. In particular, it can be shown (viz. [2] Prop. 7.1) that if a and ctq are continuous in 
Q, then one has 

/ n \i/2 

h{x, t;") = a(x) j ^ (T'Jvy] j , \Du\ - a.e. in Q (21) 

for every Borel set Q and u G BV{Q). 
Following [1] and [5], we let 

X -.= {3 e L°^{Q,R'') -.div Be L'^iQ)}. 

As proven in [5], Theorem 1.2, if denotes the outer unit normal vector to dQ, then for 
every B E X there exists a unique function [B ■ z/j^] G -Z^^n-i (c^^) such that 

f [B ■ vn]udK'-^ = [ uV-Bdfx+ [ B-Vudfx, \fueC\n). (22) 
Jdn Jn Jn 

Moreover, for u G BV{Q) each such B G L°°(f2,]R") with V ■ -B G L"(fi) gives rise to a 
Radon measure on Q, denoted {B ■ Du), satisfying the following 



[ [B ■ vn]udW-^ = f uV ■ Bdfi + [ {B ■ Du), Wu e BV{n), 
Jdn Jn Jn 



(23) 



We refer the interested reader to [H [5] for a proof. 

We shall need the following lemma, a proof of which follows from (!23|) . and the fact that 

BVioci^ \s)r\ L°^(fi) c BV{n) 

which can be easily verified. 

Lemma 3.1. Let S be as defined in Ip^) . Then 

[ [B ■ vn]udW^^ = f uV ■ Bdfi + [{B- Du) (24) 
Jdn Jn Jn 

for all u G BVioc{^\S) f] L°°(^]). 
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We conclude with a useful co-area formula for functions of bounded weighted total vari- 
ation. Details can be found in P|. 

Theorem 3.2 (Generalized Co-Area Formula). Letu G BV{Q) and suppose ^{^^^{^(^{u = 
s}) < oo holds for all s G M. Let PfpiA, Q) denote the perimeter of the set A G Q given by 



P^{A,n) := / \Dxa\ 
Jq 



Then 



/ \Du\^= / P^{{u> s},n)ds 
Jn Jr 

We note that this may, on using the representation formula f|T9l) . be recast as 



Du\^= / / h{x,v')dW"\x)ds (25) 



where is a unit outer-oriented normal vector to 1] fl d*{u{x) > s}. 



3.2 Modeling Regions with Zero or Infinite Conductivity 

Here we discuss how to formulate a suitable version of the conductivity equation ([5]) in the 
presence of inclusions of infinite and/or zero conductivity. 

Let Ooo be an open subset of fl satisfying Ooo C Q, meant to model perfectly conducting 
inclusions, and Oq be an open subset of Q with Oq C Q, meant to model insulating inclusions. 
We assume Ooo H Oq = Q \ Ooo U Oq is connected, and the boundaries 90 oo, dOo are 
piecewise C^'" for a > 0. We also assume that Oq is a mutually disjoint union of finitely many 
C^— diffeomorphic images of the unit ball, possibly empty. In addition, in two dimensions 
we require that Oq has at most one such component. 

Let a^'' and a^^ be symmetric positive definite matrix functions. For k > consider the 
conductivity problem 



m 



= on dO, 

Uk\dn = f' 



n\0o 



05 



(26) 



The perfectly conducting inclusions occur in the limiting case A; — )■ oo. The limiting 
solution is the unique solution to the problem: 



Vu = 0, 



< 



u\ 



u\ 



in fi\OoUOoo 
in Ooo 

on d{OoUOoo) 



Iao,.^i:\+dS = 0, n = l,2,... 



0, 

. u\dn = /, 



(27) 



on dOr, 
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(see the Appendix for more details), where O^o = UJ^^OJ^ is a partition of Oqo into connected 
components. Here , as in the rest of the paper, u is the outward unit normal vector and the 
subscripts ± indicate the limit taken from the outside and inside the domain, respectively. 

Remark For Lipschitz continuous conductivities in any dimension > 2, or for essentially 
bounded conductivities in two dimensions, the solutions of the conductivity equation satisfy 
the unique continuation property (see, |9] and references therein). Consequently the insulated 
(and possibly perfectly conducting) inclusions are the only open sets on which the interior 
data l^lo--! niay vanish identically. However, in three dimensions or higher it is possible to 
have a Holder continuous a and boundary data / that yield u = constant in a proper open 
subset Os C f2, see |36l|26]. We call such regions Os singular inclusions. On the other hand, 
we will not use Ohm's law in the classical sense inside perfect conductors: the current J 
inside perfectly conducting inclusions is not necessarily zero whereas Vm = within such 
regions (see |i| 122]). 

4 Uniqueness and Determining the Conformal Class 

From now on we assume that a G C"(i7, Mat(R, n)) for a > and satisfies 



where c(x) G L'^{Q \ (Ooo U Oq)) is a real, scalar- valued function, bounded away from zero 
and finite on i7\ (OooUOq) and where ctq ^ C^{^, Mat{M.,n)) is symmetric, positive definite. 

We will prove that the shape and locations of the perfectly conducting and insulating 
inclusions and the conductivity a outside of the inclusions are determined from knowledge 
of the boundary voltage /, cxo and 



a = \J(Jq^J ■ J = \ J\^-i in 

where J is the current density vector field generated by imposing the voltage / at dVL. To 
formulate our results, we first need to extend the notion of admissibility. 

Definition 4.1. A triplet of functions (/,ao,a) G H^/^{d9) x C^{n,Mat{^,n)) x L°°(^]) is 
called admissible if there exists a matrix valued function a satisfying (EHP and a divergence 
free vector field J G (L°°(fi))" such that the following three statements hold. 




(28) 




a 



|j| -1 m n. 



(29) 



(a) The vector field J satisfies 




in ^]\(OooUOo). 
in Oq, 



where u is the corresponding solution of (EZP- 

(Hi) The set of zeros of a can he decomposed as follows 




(30) 
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where Os is an open set (possibly empty), T is a Lebesgue-negligible set, and T has empty 
interior. 

We are now ready to state our main uniqueness results. 

Theorem 4.2. Let Q C M", n > 2, be a domain with connected Lipschitz boundary and let 
(/, (7o,a) G H^^'^{dQ) X C^{Q, Mat{R,n)) x L°°{Q) be an admissible triplet generated by an 
unknown conductivity cr. Define (t>{x,C,) = a{x)\C,\ao on Q\S. 

Then the potential u, solving ( 1271) . is a minimizer of the problem 

u = argmin{ / \Dv\(p : v E BV{il \ S) and vIqq = /}, (31) 
Jn 

and, if u is another minimizer of the above problem, then u = u in Vt\{\J\ = 0}. 
Moreover the zero-sets of a and |Vm| can be decomposed as follows 

{x e Q, a{x) = 0} u {x e Q : Vu = 0} =: z u r, 

where T has measure zero and Z = O^o U Oq U is open. 
Consequently 

a = — ^aoGC"(^]\Z) 

is the unique conductivity generating the admissible data triplet (/, (Jo,a). 

Remark The above theorem allows us to identify the potential u and the conductivity a 
outside the open set Z = O^o U Oq U S*. To determine if an open connected component O 
of Z is a perfectly conducting inclusion, an insulating inclusion, or a singular inclusion we 
proceed as follows: 

• If Vm = in O and a{x) ^ for some x G O, then O is a perfectly conducting 
inclusion. 

• If a = in O and u ^ constant on 90, then O is an insulating inclusion. 

• If a = in O, M = constant on dO, and a is not C° at x for some x G O, then O is 
either an insulating inclusion or a perfectly conducting inclusion. 

• If a = 0, M = constant on 50, and a G C*°((90), then the knowledge of the magnitude 
of the current {f,cro, a) is not enough to determine the type of the inclusion O. 



Proof. Suppose that u G BViod^ \ First note that for every x G \ (Ooo U Oq) there 
exists e > such that B{x, 2e) C Q and 



B{x,e) JB{x,e) 
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where J is the current density vector field described in definition ( 14. ip . Therefore 



h{x,v^)>-J-v^, \Du\-a.e. in ^]\(OooUOo). 

Thus, on using Lemma [3TT] and the fact that the current density is divergence- free away from 
the inclusions we have 



/ 

Jn 



\Du\^ = I h{x,v'')\Du\ 



n\{OouOoo) Jn\(OouOoo) 

>- [ J-v^\Du\ 

J-Du 

n\{OoUOaa) 

n-1 



J ■ vfdH' 



dn\{Oouo 
\Du\ 

Q\(OoUOoo) 



Hence 



h{x,v^) = - J ■v'', \Du\-a.e. in \ (Oq U Ooo) (32) 



Then again, since a is continuous in \ (Oq U Oqo) equation fl^Tl) gives 



But then, on f2 \ (Oq U Oqo), we have 

, 1/2 

hix^v"^) = aix)iJ2<'^>l' 

= c(x)|VM|^Jt;"|^o 

> c{x)\{Vu,v%,\ 

> aVu ■ 
= -J-v"". 

Thus it follows from (l32l) that 

J Vm 



w", \Du\-a.e. in fi\(OoUOoo). 



|J| |Vm| 

It follows from an argument similar to that of Theorem 3.5 in [28] that u = u a.e. in Vt. □ 
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5 Geometrical Properties of Equipotential Sets 

In this section we prove the area- minimizing property of the equipotential sets u~^{X) := 
fl n {x; u{x) = A} for solutions u{x) of the equation (!27|) . This generalizes results in [33] and 
[29] . The basic argument goes back to ([lOj). We shall need the following lemmas, which 
require the generalized coarea formula fl25p . 
Define the functional J-" by 

T[u] := / \Dv\^ 
Jn 

for the positive definite, matrix-valued ctq £ C^'"(fi, Mat(]R, n)). As well, for a given function 
w on Q, we introduce the notation 

w\^e '■= - min{e, max w — X,0} 

for fixed A and e > 0. For functions w on dfl, w\^^ denotes the trace at the boundary of such 
functions. 

Lemma 5.1. Let a E LipiVt), fix \ eM. and let e > 0. Then if 

u E argmin{ / |-Df I,/,; v E W]^^{yL)} 

■"\dn=f 

Then also 

MA,e E argmin{ / |-Df |<^; v E Wl_'^{Q)} 
where u\^^ = ^ min{e, max -u — A, 0} . 

Proof. For u E Wl_'^{Q) we have t>" = /i — a.e. Set u+ = maxj-u — A, 0} and = u — «+. 
Then since a and o"o are continuous we have, with 0(0;,^) = a|^|(jo) that 

j \Du\^ = [ a{x){aov'' ■ v'')\Vu\ = J^[u] 
Jn Jn 

Then again, by ( l25l) . we have 

J^[u]= [ [ a{x)d'H''-\x)dt + [ [ a{x)dK'~\x)dt 

Jt<\Ju-^{t) Jt>\Ju-^{t) 

= [ { [ a{x)Vr~^{x) + f a{x)dK'-\x)}dt 

JteR Ju~^(t) JuZ^(t) 

= + J^[n_] (33) 
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Suppose that u G argmin{J^[v];v G whereby one has with a test function w G 

Lip{Q) chosen such that w \qq = that 

T[u+] = J^[u] - J^[m_] 

< T[u + w]- J^[m_] 

< + u_ — n_ + w] 

so that M+ G ar5fmm{J^[t>]; t> G 11/^'"'^ (fi)} and, hkewise m_ G ar5fmm{J^[t>]; t> G 

^len=/+ '"\an=f~ 
follows mutatis mutandis. The lemma is immediate from two applications of this fact. □ 

Lemma 5.2. Let a G Lip{Q) and u G Lip{Q) be such that 

{x; \Vu\ = 0} = ZUL 

where Z is open, a{Z) = 0, fi{L) = and 

aVu 



G w''\n/z) 



Then 



limJ^[nA,e] = / a{x)dn''~'ix) (34) 



holds for almost every A G M. 

Proof. The proof is similar to that of Lemma 4.3 in [29]. □ 
We consider the functional 

A{T.) = jjJl-.dS (35) 

for E C n and dS the induced Euclidean surface measure on the n — 1 dimensional subset S. 
We remark that when f2 is equipped with the data dependent metric = (|o"o| |^P-i) "^(o"^^)ij 

discussed in Proposition 12.41 that the invariant volume form ^/\g\dx^ A ■ ■ ■ Adx^ on f2 induces 
the invariant volume form \ J\^-idS{x) on smooth n — 1 dimensional hypersurfaces S. For 
this reason we refer to f l35|) as an area functional. 

We are ready to establish the main result of this section, which says that equipotential 
hypersurfaces of solutions to f l27|) are minimizers of the area functional fl5S]) . 

Theorem 5.3. Let Q C M", n > 2, be a domain with connected Lipschitz boundary and 
let (/, cro,a) G C^'^(9n) x C'^'^{Q, Mat{R,n)) x L^(fi) be an admissible triplet generated by 
an unknown conductivity o G C^''^(i7 \ (Ooo U Oq) with b G (0, 1). Let v G C'^'^iVt) satisfy 

^ Ian = / and _ 

{x]\Vv\ = Q} = Z^n L^, a(Z„) = 

for Z-u open and /i(L^) = 0. 
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Then, when u is the solution to the BVP (21) we have 



A{u-\X)) < A{v-\\)) (36) 

for almost all A G M. 

Proof. Notice that when A ^ Range{u), the result is immediate. Otherwise since, under 
these hypotheses on the conductivity, u satisfies the maximum principle and, by assumption 
u and V agree on dQ, Range{u) C Range{v). If A G Range{u) / {u{Z) Uv{Z^)} is arbitrary, 
u{Z) U v{Zy) is at most countable, and since Vm and Vf are both non-zero away from 
Zand Zy respectively, u~^{\) and f ~^(A) are C^— smooth oriented hypersurfaces whereon the 
Hausdorff measure coincides with the standard Lebesgue measure. Thus 'H"'(m~^(A)) = for 
a.e. A and, by assumption 'H"~^(u^^(A) fl dVL) = for a.e. A since u extends continuously to 
the boundary, which has finite Hausdorff measure. 

Since u and v agree on the boundary ux^e \qq = vx,e \qq as well. But then by Theorem 
14.21 we have 

The result then follows from flM|) . 

□ 



6 Appendix: Perfectly conductive and insulating in- 
clusions 

The results presented in this appendix give a precise definition, by a limiting procedure, of 
potentials corresponding to conductivities that can vanish or be infinite in certain regions. 
They are slight generalization of the ones in |7] to include both perfectly conductive and 
insulating inclusions. 

Let Ooo = W^^Oi^ be an open subset of fl with O^c C Q to model the union of the 
connected components O^, (j = 1,2, ...) of perfectly conductive inclusions, and let Oq be an 
open subset of fl with Oq C fl to model the union of all connected insulating inclusions. Let 
XOoo Xoo be their corresponding characteristic function. We assume that Ooo fl Oq = 0, 
fl \ Ooo U Oq is connected, and that the boundaries dO^o, dOo are piecewise C^'" for a > 
0. Let (Ti G C"(Ooo,Mat(M,n)), and a G C"{fl \ Ooo U Oq, Mat{R, n)) be matrix-valued 
functions such that 

for constants < m, M < oo and all ^ G M". 

For each < A; < 1 consider the conductivity problem 

1 Ou 
V ■ ixuijf^i + cr)Vu = 0, — = on dOo, and uIqu = f. (37) 
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The condition on OOq ensures that Oq is insulating. It is well known that the problem ([3 
has a unique solution Uk G H^{Q) which also solves 



( V ■ aVuk = 0, 
V • cTiVwfc = 0, 

UkL = Uk\-, 



in(]\OooUOo, 
in Ooo, 
on dOr^, 



au \+ ' 

. Uk\dn = /• 



Moreover, the energy functional 



(iaiVwfc) ■ i^\_= (o-Vwfc) ■ iy\^ , on 90, 



OO ; 

on dOo, 



(38) 



h[v] = — 



Oo 



\Vv\i^dx + - 



(39) 



n\uuv 



has a unique minimizer over the maps in H^{Q) with trace / at dQ which is the unique 
solution Uk of fl38|) . 

We shall show below that the limiting solution (with A; — > 0) solves 



infi\OooUOo, 
in Ooo, 
on (90«3, 



(40) 



' V ■ crVuo = 0, 
Vuo = 0, 

Uo\+ = Uo\-, 

lao^ (^Vm) ■ iy\+ds = 0, j = 1, 2, 
tfU = 0, onSOo, 

By elliptic regularity Uq G C*^'"(fi\Ooo UOq) and for any O^'" boundary portion T of d{0oo U 
Oo), noGC'i'-((fi\(OooUOo))UT). 

Proposition 6.1. The problem [JU\ ) has a unique solution in H^{Q). This solution is the 
unique minimizer of the functional 



(41) 



'n\OooUOo 

over the set Aq := {u G H^{Q \ Oq); u\Qn = f, Vu = inOoo}- 



Proof: Note that Aq is weakly closed in H^{Q\Oq). The functional Iq is lower semicon- 
tinuous and strictly convex and, as a consequence, has a unique minimizer Uq in Aq. 

First we show that Uq is a solution of the B VP (HUl) . Since Mq minimizes (HTl) , we have 



for all (/) G i7^(r]\0o), with 







ctVmo ■ V(f>dx, 



(42) 



n\OoUOo 



we 



get li 



0, and V<^ = in Ooo- In particular, if G -ffo(f2 \ Oq 



n\OooUOo 



(V ■ (T'VuQ)(f)dx = and thus Uq solves the conductivity equation in f j40l) . If 
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we choose G H^{Q \ Oq), with = 0, and = in Oqo, from Green's formula apphed 
to we get JqoS^Wo) ■ y\+(t> = 0, V0|aoo e H^'\dOo), or equivalently, ^|aoo = 0. If 
we choose 0j G HliVt \ Oq) with 0j = 1 in the connected component O^, of Ooo and 0j = 
in Ooo \ O^, from Green's formula applied to (112!) we obtain J^^j (ctVmo) • = 0. □ 
Next we show that the equation (HUj) has a unique solution and, consequently, Uq = 
Uo\q\q^. Assume that and are two solutions and let u = U2 — Ui, then uIqq = and 



+ /aooo^^^^) ■ ^"1+"^^ + /n\o;;uo;; |Vu|2rfx = J^^o^jjo, iV^l^dx. (44) 

Thus I Vm| = in r2 \ Oq. Since \ Oq is connected and u = at the boundary, we conclude 
uniqueness of the solution of the equations fHOj) . 



Theorem 6.2. Let Uk and Uq be the unique solution of respectively in H^{Q). 
Then Uk ^ u and, consequently, Ik[uk] > Iq[u]. 

Proof: We show first that {uk} is bounded in H^iVt) uniformly in G (0, 1). Since 
1/A; > 1, we have 



A 2 1 



\^Uk\ldx+^ I \Vuk\l^dx 



n\OooUOo 



2 II ^ "'llL2(n\Oo)' 



or 



l|V«fclli2(n\Oo)^X"^''"'^(^\«o)- ^^^^ 
From fH5|) and the fact that Wfclsn = /, we see that {uk} is uniformly bounded in H^{Q\Oo) 
and hence weakly compact. Therefore, there is a subsequence Uk u* in /7^(^7\Oo), for 
some u* with trace / at dVL. 

We will show next that u* satisfies the equations fHU]) . and therefore u* = Uq on VL \ Oq. 
By the uniqueness of solutions of (jlDl) we also conclude that the whole sequence converges 
to u. 

Since Uk u* we have that = Jq^ Qouo ^^'^k ■ Vcpdx — )■ /n\ o uVm* ■ V0(ix, for all 
G C^(fi \ Ooo U Oq). Therefore V -aVM* = in f2 \ Ooo U Oq. Further, since Uk minimizes 
I[uk\ we must have Vm* = in Ooo- To check the boundary conditions, note that, for all 
G C^iyt) with = in Ooo, we have /^^^(crVMfc) ■ i/|+0(is = 0. Using the fact that were 
arbitrary, by taking the weak limit in A; — 0, we get = on dV . A similar argument 

applied to G C'^iVt) with = in Oq, = 1 in O^,, and = in Ooo\0^, also shows that 
Jg^j (ctVm*) ■ z/|+0(is = 0. Hence u* is the unique solution of the equation fHU]) on fi\Oo. 
Thus Uk converges weakly to the solution uq of ( l40l) in i7\0o. 
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7 Conclusions 



We have considered the reconstruction of an anisotropic conductivity conformal to a known 
(To when one has knowledge of the internal functional \J a^^J ■ J. Such data can be obtained 
by a novel combination of Current Density and Diffusion Tensor measurements. We have 
identified a variational problem defined in terms of the measured data and shown how to 
calculate the conformal factor from its unique solution. Further, we have presented a solution 
of the problem which allows for regions of infinite or zero conductivity, and which does not 
explicitly use Ohm's law in such regions. We also proved that the equipotential sets minimize 
the area functional corresponding to a Riemannian metric defined from the measured data. 
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